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1 INTRODUCTION 

The Surface Radiation Budget (SRB) project is a 
significant component of the Global Energy and 
Water Cycle Experiment (GEWEX) ( e.g ., Pinker et 
at., 2003; Randall et at., 2003). Targeted at its ob- 
jective, SRB has produced and archived a global 
21.5 year (July 1983-December 2004) climatology 
of, among a long list of other related variables, 
shortwave/longwave (SW/LW) downward/upward 
radiation fluxes at the Earth’s surface using the al- 
gorithm of the NASA World Climate Research Pro- 
gramme (WCRP) (Versions 2.5 and 2.6) which in- 
clude the Global Earth Observing System (GEOS) 
Version 4.0.3 meteorological information in its in- 
puts (Suarez, 2005). (The earlier versions of these 
data sets have been studied from different per- 
spectives. (e.g., Gupta et at., 1999; Zhang et at., 
2003; Zhang et at., 2004; Pinker, Zhang and Dut- 
ton, 2005). 

The uncertainties in the atmospheric information 
and imperfect algorithms in radiation transfer mod- 
els make the satellite-derived surface information 
uncertain. Therefore, ground-based observations 
are indispensable for corroboration and ultimate 
validation of the surface information derived from 
satellite-observations. At the same time, surface 
observations can be used to train radiation transfer 
models. Additionally, the surface observations are 
more reliable for detecting long-term trend related 
to climate change. 

For this reason, the GEWEX program has a 
Baseline Surface Radiation Network (BSRN) de- 
signed to specialize in ground-based observations. 
In this study we use 2849 site-months of BSRN data 
representing 35 ground stations. 
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In addition, we use the World Radiation Data 
Centre (WRDC) and Global Energy Balance 
Archive (GEBA) data. The WRDC data used here 
include daily means of shortwave downward fluxes 
from 474 sites of the World Radiation Data Cen- 
tre (WRDC) depository. The GEBA data include 
monthly means of shortwave downward fluxes from 
797 sites and longwave from 20 sites. 

In this investigation, we make systematic SRB- 
BSRN, SRB-WRDC and SRB-GEBA comparisons 
for both shortwave and longwave daily and monthly 
mean radiation fluxes at the Earth’s surface. 

We first have an overview of all the compara- 
ble pairs of data in scatter or scatter density plots. 
Then we show the time series of the SRB data at 
grids in which there are ground sites where long- 
term records of data are available for comparison. 
An overall very good agreement between the SRB 
data and ground observations is found. 

To see the variablity of the SRB data during the 
21 .5 years, we computed the global mean and its 
linear trend. No appreciable trend is detected at the 
5% level. 

The empirical orthogonal functions (EOF) of the 
SRB deseasonalized shortwave downward flux are 
computed over the Pacific region, and the first EOF 
coefficient is found to be correlated with the ENSO 
Index at a high value of coefficient of 0.7083. 


2 Data 


The satellite data is the SRB data, and the ground- 
based data are from the BSRN, the WRDC and the 
GEBA databases. 
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Northern Latitude (°) 


Figure 1: Scaled grid areas in each 1 “-latitudinal 
belt. The globe is first divided into 1 80 1 “-latitudinal 
belts. In Northern Hemisphere, from 0°N to 45°N, 
each belt has 360 grids; from 45°N to 70°N, each 
belt has 1 80 grids; from 70°N to 80°N, each belt has 
90 grids; from 80°N to 89°N, each belt has 45 grids; 
from 89°N to 90°N, the belt has 3 grids. The South- 
ern Hemisphere is a mirror image of the Northern 
Hemisphere about the Equator. 

2.1 The SRB Data 

The SRB data are output of models designed for 
the SRB project, and the data sets are presented 
for a network of 44,016 grids of nearly equal area 
as shown in Figure 1. The Earth’s surface is first 
divided into 1° latitudinal belts, and each such belt 
is then divided into grids of equal area, the number 
of grids depending on the latitude. In the tropical 
zone, the grids are 1° x 1° and, the northernmost 
and southernmost belts are divided into 1° x 120° 
grids. 

2.2 The BSRN Data 

We use the BSRN archived data worth 2849 site- 
months, representing 35 ground stations ranging 
from the northernmost one at Ny Alesund in Spits- 
bergen (78.93° N, 11.95°E) to the southernmost 
one at South Pole (90° S), covering the time pe- 
riod from January 1992 to August 2004. The ob- 
served variables used in the current study are short- 
wave and longwave downward fluxes at the Earth’s 
surface. The BSRN observations are made at 1- 
, 2, 3- or 5-minute intervals. More information 
about BSRN can be found in Ohmura et al. (1998); 
Philipona et al. (1998); McArthur, 1998; Wild et 
al. (2005). Figure 2 shows the locations of the 35 
BSRN sites. 



Figure 2: Baseline Surface Radiation Network 

(BSRN) sites with data. 


The data are first processed to generate 3-hourly, 
3-hourly-monthly, daily and monthly means so that 
they can be compared with their BSRN counter- 
parts on the same time scales. Note that the SRB 
3-hourly data represent instantaneous snapshots of 
the Earth at UTC 00:00, 03:00, 06:00, ..., 21:00 
of averages over grids of about 11,588 Km 2 , and 
these hours are called “time stamps”. The BSRN 
3-hourly means, on the other hand, are averages of 
all valid values observed at 1-, 2-, 3- or 5-minute in- 
tervals during 3-hour periods centered on the said 
“time stamps”, and these averages represent sin- 
gle points on the ground. Contrary to the intuition, 
this approach actually makes the SRB 3-hourly data 
and corresponding BSRN 3-hourly means more 
comparable, as shall be made clearer later. 

2.3 The WRDC Data 

The World Radiation Data Centre (WRDC), estab- 
lished in accordance with Resolution 31 of WMO 
Executive Committee XVIII in 1964, is a central 
depository of solar radiation data collected over 
1000 ground sites over the world. In this study, 
we use a subset of the WRDC data which in- 
cludes daily mean shortwave downward flux from 
474 sites throughout the world spanning the years 
1983-1993. Out of these sites, 55 have continuous 
records, and the records at other sites are sporadic 
at various extent. Monthly means are computed 
whenever it is possible. Figure 3 shows the 474 
WRDC sites. 

2.4 The GEBA Data 

The Global Energy Balance Archive (GEBA) is 
Project A7, Water (WMO/ICSU), of the World Cli- 
mate Programme. In this study, we use the monthly 
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World Radiation Data Center (WRDC) Sites with Any Daily Mean 
Shortwave Downward Flux Data from 1983 to 1993 


Longitude (°) 


Elevation (m) 

Figure 3: World Radiation Data Center (WRDC) 
sites with any data from 1 983 to 1 993. 


Figure 4: 797 Global Energy Balance Archive sites 
with data from 1 983 to 2003. 

means of shortwave downward fluxes from 797 
GEBA sites and longwave from 20 GEBA sites. Fig- 
ure 4 shows the shortwave GEBA sites. 


3 The Comparison 

Before we look at those ground sites with long-time 
records, we first have an overview of all available 
BSRN data points paired with their SRB counter- 
parts in plots of scatters or scatter density, which 
include all sporadic yet valid data points. 

3.1 An Overall View 

Figure 5 shows the scatter density plot of daily 
mean SRB Version 2.6 shortwave downward fluxes 
versus their BSRN counterparts for the years 1 992- 
2004. The 35 BSRN sites represented are shown 


in Figure 2. The reason that a scatter density plot 
is chosen over a density plot is that significant over- 
lapping occurs in a plot of the latter kind due to the 
large number of data points. The color of a grid in- 
dicates the relative count of the data points in the 
grid which is the absolute count of the data points 
divided by the total, i.e., N. 

Figure 6 is the same as Figure 5 except it is the 
monthly mean and scatters in place of scatter den- 
sity are plotted. 

Figures 7 and 8 are the longwave counterparts of 
Figures 5 and 6. 

Figures 9 and 10 are the same as Figures 5 and 
6 except the ground-based observations are from 
the WRDC database. Note the large number of 
data points, N = 1640247 , representing 655 differ- 
ent WRDC ground sites from 1983-1993 (Figure 3, 
in the daily means and the concentration of high 
values of density along the white line, the ideal sit- 
uation. The derived monthly means, representing 
474 ground-sites, still have a number of data points 
large enough to justify a scatter density plots. 

The GEBA data are available only in monthly 
means. Figure 11 shows the monthly mean GEBA 
shortwave downward fluxes versus their SRB coun- 
terparts for the years 1 983-2003. Both the numbers 
of GEBA ground sites and data points are larger 
than that of WRDC data. The 82,977 data points 
represent 797 ground sites as shown in Figure 4. 

Only 19 of the GEBA sites have 1013 available 
monthly mean longwave downward flux data points 
which are shown in Figure 12 along with their SRB 
counterparts. 

Generally speaking, the updated satellite-derived 
SRB data are in very good agreement with ground 
observations. In terms of the Bias absolute value 
and RMS, the longwave data show much better 
agreement than the shortwave data for understand- 
able reasons, and of the shortwave data, the GEBA 
data, with the largest numbers of ground sites and 
data points of all the three sets of ground-based 
data, show the best agreement with the SRB data. 

3.2 Time-Series of Monthly Means 

Fortunately, all the three sets of ground-based data 
have sites with considerably long records of obser- 
vations. 

Figures 13, 14 and 15 show the time series of 
monthly means of shortwave downward fluxes at 
BSRN sites NYA (Ny Alesund, Spitzbergen), KWA 
(Kwajelein, Marshall Islands) and SPO (South Pole) 
for the years 1 992-2003 along with their SRB coun- 
terparts. 



797 Global 


Sites with Data from 1 983-2003 


Longitude (°) 
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Scatter Density of Daily Means of SRB(V2.6)-BSRN SWDW Fluxes 
for 1992 - 2004 



BSRN SWDW (W m' 2 ) 


Figure 5: Scatter density of SRB(V2.6)-BSRN 

shortwave downward flux daily comparison.The 
white line represents the points where BSRN val- 
ues are equal to their SRB counterparts, and the 
black line is the result of linear regression. R 2 is 
the coefficient of multiple determination; Bias is the 
average of the SRB-BSRN difference, and RMS is 
the root mean square of the same difference; p is 
the SRB-BSRN correlation coefficient; a is the stan- 
dard deviation of the SRB-BSRN difference; pbsrn 
is the average of all BSRN data points; N is the total 
number of data points in the plot. The same nomen- 
clature is used throughout the paper. 

Comparison of SRB(V2.6) and BSRN Data for All BSRN Sites for 1992 - 2004 



Figure 6: SRB(V2.6)-BSRN shortwave downward 
flux monthly mean comparisons. 


Scatter Density of Daily Means of SRB-BSRN LWDW Fluxes 
for 1992 - 2001 



Figure 7: Scatter density of SRB(V2.5)-BSRN long- 
wave downward flux daily mean comparison. 
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Figure 8: SRB(V2.5)-BSRN longwave downward 
flux monthly mean comparison. 
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Scatter Density of Daily Means of SRB-WRDC SWDW Fluxes 
for 1983 - 1993 
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Figure 9: Scatter density of SRB(V2.5)-WRDC 
shortwave downward flux daily mean comparison. 
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Scatter Density of Monthly Means of SRB-WRDC SWDW Fluxes 
for 1983 - 1993 



Figure 10: Scatter density of SRB(V2.5)-WRDC 
shortwave downward flux monthly mean compari- 
son. 




Figure 13: Time series of monthly means of 

SRB(V2.6) and BSRN shortwave downward fluxes 
and difference thereof at BSRN site NYA (Ny 
Alesund, Spitzbergen; 78.93°N, 1 1 .95° E). 


Scatter Density of Monthly Means of SRB(V2.6)-GEBA SWDW Fluxes 
for 1983 - 2003 



Figure 11: Scatter density of SRB(V2.6)-GEBA 
shortwave downward flux monthly mean compari- 
son. 

Scatter Density of Monthly Means of SRB(V2.5)-GEBA LWDW Fluxes 
for 1983 - 2003 



Figure 12: Scatter density of SRB(V2.5)-GEBA 
longwave downward flux monthly mean compari- 
son. 



Year 


Figure 14: Time series of monthly means of 

SRB(V2.6) and BSRN shortwave downward fluxes 
and difference thereof at BSRN site KWA (Kwa- 
jalein, Marshall Islands; 8.72°N, 167.73°E). 


Figures 16, 17 and 18 show the longwave coun- 
terparts of the above three BSRN sites. 

Eight of the 35 BSRN sites have continuous 
records at the same time span. 

Figures 19, 20, 21 and 22 show the time se- 
ries of WRDC monthly mean shortwave downward 
fluxes for the years 1 983-1 993 at WRDC sites 0626 
(Wrangel Island, Russia), 0350 (Klote, Switzer- 
land), 0857 (Narok, Kenya), and 0953 (Bulawayo, 
Zimbabwe) along with their SRB counterparts. 

Out of the 655 WRDC sites represented in this 
study, 55 have continuous records over the same 
time span. 

Figures 23, 24, 25 and 26 show the time series 
of the GEBA monthly mean shortwave downward 
fluxes for the years 1983-2003 at GEBA sites 1412 
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Figure 15: Time series of monthly means of 

SRB(V2.6) and BSRN shortwave downward fluxes 
and difference thereof at BSRN site SPO (South 
Pole; 90.00°S, 0.00°E). 



Year 


Figure 18: Time series of monthly means of 

SRB(V2.5) and BSRN longwave downward fluxes 
and difference thereof at BSRN site SPO (South 
Pole; 90.00°S, 0.00°E). 



Figure 16: Time series of monthly means of 

SRB(V2.5) and BSRN longwave downward fluxes 
and difference thereof at BSRN site NYA (Ny 
Alesund, Spitzbergen; 78.93°N, 11.95°E). 




Figure 19: Time series of monthly means of 

SRB(V2.5) and WRDC shortwave downward fluxes 
and difference thereof at WRDC site 0626 (Wrangel 
Island, Russia; 70.58°N, 178.32°W). 
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Figure 17: Time series of monthly means of 

SRB(V2.5) and BSRN longwave downward fluxes 
and difference thereof at BSRN site KWA (Kwa- 
jalein, Marshall Islands; 8.72°N, 167.73°E). 





SRB(V2.5) and WRDC shortwave downward fluxes 
and difference thereof at WRDC site 0350 (Klote, 
Switzerland; 47.29°N, 8.32°E). 
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Year 


Figure 21 : Time series of monthly means of 

SRB(V2.5) and WRDC shortwave downward fluxes 
and difference thereof at WRDC site 0857 (Narok, 
Kenya; 1.08°S, 35.50°E). 



Figure 22: Time series of monthly means of 

SRB(V2.5) and WRDC shortwave downward fluxes 
and difference thereof at WRDC site 0953 (Bul- 
awayo, Zimbabwe; 20.09°S, 28.37°E). 




SRB(V2.6) and GEBA shortwave downward fluxes 
and difference thereof at GEBA site 1412 (Kiruna, 
Sweden; 67.85°S, 20.23°E). 


(Kiruna, Sweden), 1238 (Jokioinin, Finland), 1176 
(Uccle, Belgium) and 1274 (Ajaccio, France) along 
with their SRB counterparts. 

Altogether, 101 of the 797 GEBA sites have 
nearly continuous records during the period of two 
decades from 1983-2003. 

The results presented in Figures 13-26 show the 
the SRB Version 2.5 and 2.6 data not only catch the 
magnitudes of the shortwave and longwave down- 
ward fluxes consistently with considerable accu- 
racy, but the seasonality as well. 

4 The Variability of the Data 

Here we examine if there is any appreciable linear 
trend in the SRB data over the course of 21 .5 years 
and how the data variability is associated with the 
El Nino Southern Oscillation (ENSO). 

4.1 Linear Trend 

In Figure 27, the green curve in the upper panel 
shows the sequence of the global mean of the SRB 
(Version 2.5) monthly means of shortwave down- 
ward fluxes at the Earth’s surface, and the red curve 
shows the average of the SRB data at selected 
grids corresponding to 101 GEBA sites with con- 
tinuous records. The purpose of the red curve is to 
see to what an extent the average of ground obser- 
vations from these sites may represent the global 
average. 

Interestingly, the seasonality of the SRB global 
mean is synchronous with that of the Southern 
Hemisphere, and the average from the selected 
grids shows a Northern Hemisphere seasonality 
which is much more pronounced than that of the 
global mean, apparently due to the fact that most of 
the 101 GEBA sites are in Northern Hemisphere. 

The lower panel of Figure 27 shows the desea- 
sonalized version of the upper panel. At the same 
time, the trends as derived from linear regression 
are shown in straight lines. Following the method 
of Weatherhead etal. (1998), autocorrelation of the 
data is considered. Computation of the partial auto- 
correlation coefficients indicates that only lag-1 au- 
tocorrelation is statistically significant, and conse- 
quently, an autoregressive model of order 1, it i.e., 
AR(1) model can be justifiably used to estimate the 
linear trend in the mean averaged both globally and 
at selected GEBA sites. 

Thus the trends computed are b S ELECTED = 
0.0454 (Wm _2 a _1 ) and bcLOBAL = 0.0447 

(Wm _2 a^ 1 ). Both the 95% confidence intervals 
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Figure 24: Time series of monthly means 

of SRB(V2.6) and GEBA shortwave downward 
fluxes and difference thereof at GEBA site 1238 
(Jokioinen, Finland; 60.82°S, 23.50°E). 
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Figure 25: Time series of monthly means of 

SRB(V2.6) and GEBA shortwave downward fluxes 
and difference thereof at GEBA site 1176 (Uccle, 
Belgium; 50.80°S, 4.35°E). 
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Figure 26: Time series of monthly means of 

SRB(V2.6) and GEBA shortwave downward fluxes 
and difference thereof at GEBA site 1274 (Ajaccio, 
France; 41 .92°S, 8.80°E). 



Figure 27: SRB(V2.5) shortwave monthly means 
averaged globally and at 101 GEBA composite 
sites. The upper panel shows the time sequences, 
and the lower panel shows the corresponding de- 
seasonalized versions as well as the linear trends 
displayed in the data. 


are wide due to the autocorrelation and stride 0, 
leaving small positive trends in uncertainty. 

4.2 Correlation with the ENSO 

To see the variablity of the SRB data, we computed 
the empirical orthogonal functions (EOFs) of the 
shortwave downward fluxes, and to emphasize the 
component associated with the ENSO, the range of 
interest is limited the Pacific region. The first EOF 
of the deseasonalized version of the data explain 
20.60% of the total variance, and its coefficient is 
shown in Figure 28 along with the Southern Oscil- 
lation Index (SOI), and the correlation coefficient of 
the two curves is as high as 0.7083. Note that the 
boundary of the Pacific region is indicated in the fig- 
ure. 

As the figure indicates, there is a significant de- 
terministic component of the shortwave downward 
flux that is associated with the ENSO. If a trend is to 
be computed for the reference of practical decision- 
makers, this process, along with other possible de- 
terministic oscillatory processes should be identi- 
fied and eliminated first. 


5 Discussion and Conclusions 

Despite the challenges involved in the GEWEX 
SRB project, the SRB Versions 2.5 and 2.6 al- 
gorithms generated the most comprehensive sur- 
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Figure 28: The coefficient of the first EOF of the 
SRB(V2.5) deseasonalized downward shortwave 
flux from July 1983 to December 2004 over the Pa- 
cific in comparison with the Southern Oscillation In- 
dex (SOI) over the same period. Note that the cor- 
relation coeffcient, p, is as high as 0.7083. 


face radiation database from satellite-based obser- 
vations. The improvement of Version 2.6 over Ver- 
sion 2.5 is that more accurate astronomical geom- 
etry of the sun is introduced. The time span of 
the database covers the time from July 1983 to 
December 2004, and the resolution is based on a 
nearly equal-area scheme of 44,01 6 grids which are 
1° x 1° around the Equator. The shortwave and 
longwave fluxes at the Earth’s surface representing 
daily and monthly means turned out to be in very 
good agreement with the independently observed 
ground-based BSRN, WRDC and GEBA data. The 
SRB 3-hourly and 3-hourly-monthly data which are 
not shown in this study are equally in good agree- 
ment with their BSRN counterparts except they are, 
understandably, noisier than the daily and monthly 
means as manifested in their larger standard devia- 
tions. 

In addition to the solar zenith angle, the cloudi- 
ness plays a significant role in the shortwave flux. 
The great temporal and spatial variablity of clouds 
manifests itself in the variablity of shortwave flux 
at the Earth’s surface. However, the noise exhib- 
ited in the scatter and scatter density plots cannot 
be equally attributed to the SRB and ground-based 
data sets. As separate computation indicates, the 
ground-based time series contains slightly more 
variance than their SRB counterparts, which is 
due to the fact that each ground-based observa- 
tion is conducted at a single point on the surface 
of the Earth that may experience the whole range 


of cloudiness from zero to the possible maximum 
while each of the SRB value represents an instant 
spatial average over a grid whose linear size ranges 
from about 111 to 192 Km. An SRB grid is thus 
less likely to experience the total cloud-free condi- 
tion. So the time averaging process has removed, 
but not completely, the extra variance of the ground- 
based data as compared to their SRB counterparts. 

Out of the ground-based sites, 8 BSRN sites, 55 
WRDC sites and 101 GEBA sites have nearly con- 
tinuous records over extended periods, and they 
are generally in good and consistent agreement 
with their SRB counterparts. In the case of SRB- 
BSRN comparisons, the agreement tends to be bet- 
ter in lower latitudes than higher ones. Although we 
tend to treat the ground-based observations as the 
“truth”, we need to be careful in evaluating the qual- 
ity of satellite-ground agreement. There may be er- 
rors attributable to technical aspects of the instru- 
ments, the micro-environments around the instru- 
ments, and human errors as well. 

Compared with the shortwave fluxes, the long- 
wave fluxes behave more predictably, exhibiting 
much less irregularity, and the scatter and scatter 
density plots show few points that can be consid- 
ered as outliers. This can be attributed to the fact 
that the longwave radiation intensity is determined 
by the temperature profile of the atmosphere and 
other components residing therein, and the temper- 
ature is not as prone to change as the solar radi- 
ation transferring in the atmosphere. As a conse- 
quence, both the bias and RMS errors in the long- 
wave comparisons are markedly smaller than their 
shortwave counterparts. 

The trend estimate on the global and regional 
scales is interesting and a challenge both mathe- 
matically and for practical decision-makers. Based 
on an AR(1) model, the global mean of the SRB 
Version 2.5 does not show an appreciable linear 
trend over the years from July 1 983 to June 2003. 

The first EOF coefficient of the deseasonalized 
global mean of monthly means of the SRB Version 
2.5 shortwave downward flux over the Pacific region 
explains a significant part of the total variance and 
is highly correlated with the ENSO. This means that 
climatic variations associated with the ENSO leave 
unmistakable signals in the time sequence of the 
global radiation field. 
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